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The metric of a perturbed Robertson- Walker spacetime is characterized by three 
functions: a scale-factor giving the expansion history and two potentials which gen- 
eralize the single potential of Newtonian gravity. The Newtonian potential induces 
peculiar velocities and, from these, the growth of matter fluctuations. Massless par- 
ticles respond equally to the Newtonian potential and to a curvature potential. The 
difference of the two potentials, called the gravitational slip, is predicted to be very 
small in general relativity but can be substantial in modified gravity theories. The 
two potentials can be measured, and gravity tested on cosmological scales, by com- 
bining weak gravitational lensing or the Integrated Sachs- Wolfe effect with galaxy 
peculiar velocities or clustering. 
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1. Introduction 

Cosmic acceleration remains as mysterious today as it was when its discovery was 
announced in 1998 (Riess et al. 1998; Perlmutter et al. 1999). In the cosmological 
standard model based on general relativity, a homogeneous, isotropic, and uniformly 
expanding universe accelerates if and only if the pressure is large and negative such 
that p + ip < 0. Dark energy is the name given to a class of exotic fluids having 
this property. Dark energy may be either constant (the cosmological constant A) or 
have varying density and pressure (dynamical dark energy). At present, the chief 
goal of observational cosmology is to measure the expansion history H{a) in order 
to determine the time-dependence of the energy density p{a) from the Friedmann 
equation 

H\a) = ^Gp{a) - Ka-K (1.1) 

Measurements of p{a) yield the dark energy equation of state p{p) through energy 
conservation, dp/dliia = —3{p + p). 

The conclusion that dark energy exists depends on the assumptions of the cos- 
mological standard model. For example, if the galaxy distribution is strongly in- 
homogeneous while preserving near-isotropy around us, it is possible for galaxy 
motions along the past light cone to mimic the acceleration of a homogenous uni- 
verse without dark energy (Kolb et al. 2005; Vanderveld et al. 2006; Clifton et 
al. 2008). This possibility requires very special initial conditions that break the 
fundamental translation-invariance of the standard model; moreover, gravitational 
instability might convert large radial gradients into large angular gradients that are 
not seen. The breaking of translation invariance would be problematic for initial 
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conditions generated by inflation, among other challenges. This article will there- 
fore retain the assumption of large-scale homogeneity and isotropy present in the 
standard cosmological model. 

There is one key element of the standard cosmological model that is poorly 
tested: general relativity (GR) itself. Since Einstein's 1917 paper introducing the 
cosmological constant, GR has been so fundamental to cosmological models that 
it is not immediately obvious how to modify it. Clearly any modifications great 
enough to eliminate the need for dark energy must preserve the successes of GR in 
the solar system and in binary pulsars (Will 2006). Thus the modifications of the 
field equations or their solution must be scale-dependent: negligible on short length 
scales and appreciable on the Hubble scale. 

Measurements of the expansion history alone cannot distinguish GR from alter- 
native theories. An example of a modified gravity theory having cosmic acceleration 
without dark energy is the DGP braneworld model of Dvali et al. (2000), imple- 
mented in cosmology by Dcffayet et al. (2002). This is a higher-dimensional theory 
in which our universe lives on a brane of three spatial dimensions. The field equa- 
tions are modified at long wavelengths. In the so-called self-accelerating branch of 
the theory, the Friedmann equation is modified to become 



where Hoc is a constant. At early times H > i?oo- As the universe expands, p{a) 
decreases so H decreases until it asymptotically approaches H^o- A constant Hubble 
parameter H = dlna/dt implies exponential expansion, the same as a cosmological 
constant-dominated universe in GR. Such a model has acceleration without dark 
energy. 

The DGP model makes a specific prediction for H{a) that can be tested and 
falsified; it appears to be in confiict with observations (Fairbairn & Goobar 2006; 
Lombriser et al. 2009). However, there are other classes of modified gravity theories 
that arc not, and cannot be, ruled out based on the expansion history. In particular, 
models in which the Einstein-Hilbert action for gravity is modified by an arbitrary 
function of the Ricci scalar R can reproduce any H{a) (Multamaki & Vilja 2006; 
Capozziello et al. 2006). The action for the theory is 



where f{R) = for GR. There is no fundamental theory for f{R); instead, there 
exist parametrized forms (e.g. Hu & Sawicki 2007a; Starobinsky 2007) for which 
predictions can be made and compared with observations. 

These results imply that GR cannot be tested on cosmological scales using the 
expansion history. One must look beyond the homogeneous and isotropic, uniformly 
expanding Robertson- Walker models to include perturbations. Moreover, it is most 
helpful to identify those features of the data that most distinguish modified gravity 
theories from GR, and use them to construct cosmological tests. Theorists and 
observers have begun to do this (e.g., Zhang et al. 2007; Reyes et al. 2010). As we 
show, these tests involve two functions of space and time, which we call potentials. 
We also show that observations of galaxies or cold dark matter, in particular the 
growth of structure, are insufficient to test GR. 
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What does it mean to test GR? General relativity makes a number of assump- 
tions that may also hold in alternative theories of gravity (Bertschinger 2006a): 

1. Spacetime is a 4-dimensional pseudo-Riemannian manifold. 

2. Special relativity holds locally. In particular, energy-momentum is locally con- 
served. 

3. The weak equivalence principle holds, i.e., freely-falling bodies follow space- 
time geodesies. 

4. The metric is the soution to the Einstein field equations subject to appropriate 
initial and boundary conditions; this is uniquely true in GR. 

The first three assumptions will be adopted here for modified gravity theories as well 
as for GR. Modifying gravity then implies modifying the field equations; the fourth 
assumption is to be tested. Specifically, testing GR requires measuring the 
metric (or the transfer functions and power spectra for perturbations in 
the metric) and comparing measurements with the solution of the field 
equations. 

Section 2 introduces the two potentials used to characterize cosmological per- 
turbations and shows that nonrelativistic matter (atoms and cold dark matter) are 
sensitive to only one of them while the deflection and redshift of photons are sen- 
sitive to both potentials. Section 3 introduces the primary tools for measuring the 
potentials: peculiar velocities, the growth of structure, the Integrated Sachs- Wolfe 
effect, and gravitational lensing. While full reconstruction of the potentials remains 
a longterm goal, current efforts emphasize measurements of correlation functions 
and power spectra, which also provide valuable constraints on the two potentials. 

Section 4 summarizes the fleld equations for the potentials in GR and in one 
class of modified gravity theories, the f{R) theories. It is shown that a key discrim- 
inating feature between theories is the "gravitational slip," the difference of the 
two potentials. In GR the slip is generated only by relativistic shear stress and is 
negligible in the matter-dominated era. In f{R) theories the slip is generated by a 
new gravitational field called the scalaron. The scalaron also modifies the growth 
of structure while avoiding solar-system constraints. Section 5 briefly summarizes 
current tests. 

More extensive reviews are given by Caldwell & Kamionkowski (2009), Silvestri 
& Trodden (2009), Uzan (2009), and Jain & Khoury (2010). 

2. Two potentials 

The most general metric of a weakly perturbed Robertson- Walker spacetime — 
regardless of whether GR is valid — takes the form 

ds^ = a2(r) {-(1 + 2^)dT^ + 2widx'dT + [(1 - 2^)jij + 2sij]dx'dx^} (2.1) 

where "fij is the three-metric of constant-curvature spaces, e.g. jij ~ Sij for Carte- 
sian coordinates in a Euclidean space. The inverse spatial metric is 7'-^ . For an 
unperturbed Robertson- Walker spacetime, $ = ^ = = Sjj = 0. Without loss of 
generality we can take to be symmetric and traceless, i.e., 7'-'sij = 0. The time 
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variable r is called conformal time and is related to the proper time for comoving 
observers (those at fixed spatial coordinates) in the unperturbed background by 
dt ~ a dr. 

Equation (|2.ip is too general in that the same spacetime is described by infinitely 
many sets of solutions for the functions ($, Wi, Sij). This ambiguity arises because 
of general covariance, i.e. the invariance of the gravitational action under coordi- 
nate transformations. Like the vector potential of electromagnetism, the metric 
perturbations change under gauge transformations (here, infinitesimal coordinate 
transformations). We are free to remove the ambiguity by gauge- fixing, which means 
choosing a particular coordinate system. We choose the Poisson gauge (Bertschinger 
1996), a coordinate system defined by transverse gauge conditions 



where spatial divergences are covariant derivatives taken with respect to the back- 
ground 3-metric 7ij. The function Wi is called the gravomagnetic potential and is 
associated with vector perturbations. The function Sij is called the gravitational 
wave strain and is associated with tensor perturbations. The Poisson gauge gener- 
alizes the transverse-traceless gauge widely used for gravitational waves to include 
all possible perturbations of the metric (scalar, vector, and tensor). 

Although gravomagnetic and gravitational wave perturbations are expected to 
be present in the universe, they have not been directly detected and are known to 
be much smaller in amplitude than the other terms in the metric. Therefore we 
neglect them in what follows. 

The potential $ is called the Newtonian potential. Being part of the time-time 
part of the metric, it is important for the motion of nonrelativistic particles. For 
slowly moving particles, dx^ /dr is small and the spatial part of the metric makes a 
negligible contribution to the proper time and therefore to the equations of motion 
that extremize proper time. The second potential 4' is called the curvature potential 
because it is associated with spatial curvature. It is important for the motion of 
relativistic particles. 

With this coordinate choice and with the neglect of gravomagnetic and gravita- 
tional wave perturbations, the metric simplifies to the Conformal Newtonian metric 
(Mukhanov et al. 1992, Ma & Bertschinger 1995) 



The potentials are assumed to be small enough ($ v]/ ^ lO"'"^) that their products 
can be neglected. For a flat background, with jij — Sij, and with 4' = $, this 
metric is conformal to (i.e., identical to aside from an overall multiplicative factor) 
the weak- field metric used in the Newtonian limit (Hartle 2003). The conformal 
Newtonian coordinate system used here has the advantage of making motion look 
simple. The reader should beware that different choices of names and signs for the 
two potentials appear in the literature. 

Freely-falling particles move along geodesies of metric (|2.3p and obey equations 
of motion 



V'u;, = , V^sy = , 



(2.2) 



ds^ = a^(T) [-(1 + 2<^>)dT^ + (1 - 2^)-f,jdx^dx^] . 



(2.3) 



dx' 

d7 



{i + ^ + ^y , 



(2.4a) 



1 d 
7a(l - ^) dr 



7a(l - '^)v^ 



V»($ -I- v'^^) - (1 $ + ^)-i]kV^v^ . (2.4b) 
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Here, is the proper (physical) 3- velocity measured by a comoving observer (note 
that r is not the proper time). Specifically, it is the peculiar velocity with respect 
to uniform Hubble fiow. Other terms in the equations of motion have their usual 
definitions: = jijV^v^ , 7 = (1 — ti^)"^/^ is the special relativistic Lorentz factor, 
and "fjj. is the spatial connection needed when non-Cartesian coordinates are used 
(giving rise to terms such as centripetal acceleration). The factors (1 + $ + ^') 
convert proper velocity to coordinate velocity. Hereafter we will ignore the ^ 10^^ 
corrections from proper to coordinate velocities. 

The equations of motion simplify in the limits u — >■ (nonrelativistic particles) 
and w — ^ 1 (massless particles). In the first case, denoted cold dark matter (CDM), 



1 d{av) 
a dr 



V$ , < 1 (CDM) . (2.5) 



In the second case. 



dv 

— = -V_l($ + *) , = 1 (photons) (2.6) 
or 

where the gradient is taken in the plane perpendicular to the photon trajectory: 

VI = V - v'vjW^ . (2.7) 

Equations (j2.5l) and (j2.6p are very important and are more general than GR. 
They have a crucial implication: The dynamics of galaxies and CDM cannot 
fully test GR because they depend on only one of the two potentials. 

Gravitational lensing measurements are sensitive to the sum of the two poten- 
tials. Another dependence arises for the redshift of light: 

d\n(aE) „ . , , , , 

= V • V$ , = 1 (photons) . 2.8 

dr OT 

The first term on the right-hand side is less familiar than the second term (gravita- 
tional redshift). It becomes more familiar when equation (j2.8l) is integrated to give 
the accumulated change in energy of a photon between emission at redshift z and 
observation: 

(1 + z)§^ = 1 + $em - $obs - / (^ + -^Kdr/dz) dz , (2.9) 

i^evci Jo 

where a dot denotes d/dr. Equation (j2.9p gives the familiar Sachs- Wolfe effect 
(Sachs & Wolfe 1967); the last term, arising from time- variation of the potentials 
along the past light cone, is called the Integrated Sachs- Wolfe (ISW) effect. Thus 
cosmic microwave background (CMB) anisotropy responds to both the Newtonian 
potential $cin and to the time derivative of the sum of the two potentials. 

3. Methods to measure the potentials 

Testing gravity on cosmological scales requires making distinct measurements for 
each of the three functions appearing in the metric. The most direct way is: 
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1. Measure a(r) using redshift and distance measurements to infer the cosmic 
expansion history H(a)\ integrate dr/da — \/{Ha?) to get r(a) which can 
then be inverted. 

2. Measure $(a;,T) using equation (j2.5p with the pecuUar velocities of galaxies. 

3. Measure (<I' + 4')(a;, r) using equation (j2.6|) with gravitational lensing or equa- 
tion (ES) with the ISW effect. 



Although the second two measurements give derivatives of the potentials rather 
than their absolute values, with appropriate initial or boundary conditions the 
potentials can be determined up to irrelevant constants. These measurements are 
then compared with the relationships between the three functions predicted by GR 
or modified gravity theories. 

In practice these direct methods are difficult at best. Peculiar velocities are dif- 
ficult to measure with sufficient precision for gravitational tests; while the radial 
velocity is easy to measure from redshifts, the Hubble velocity is proportional to 
distance. Even with the precision given by Type la supernovae, uncertainties in 
distance measurements exceed the peculiar velocities for distances beyond a few 
tens of Mpc. Gravitational lensing is promising, however the deflections are small 
for structures larger than a few Mpc. The ISW signal is concentrated at the lowest 
multipoles of the CMB where primary anistotropy and cosmic variance (the vari- 
ance arising from Gaussian random processes with a very small number of terms 
contributing to the signal) are large. 

(a) Growth of structure 

To increase the power of galaxies to measure $, cosmologists commonly turn to 
the growth of density perturbations. Density perturbations are part of the stress- 
energy tensor which is decomposed for scalar perturbations as 

T\ = -{p + 5p), (3.1a) 
= -ip + p)V,u, (3.1b) 

TV = S^^{p + Sp) + l{p + p)(^W,~l6^^Ay (3.1c) 

where A = 7'-' ViVj is the spatial Laplacian, p andp are the background density and 
pressure and u and tt are velocity and shear stress potentials, respectively. The per- 
turbations are gauge-dependent. A convenient gauge- invariant density perturbation 
is the physical number density perturbation in the fluid rest frame (Bertschinger 
2006b), obtained from conformal Newtonian gauge variables as 

V = — ; ; — , n = - = atl . (3.2) 

p + p a 

For a nonrelativistic fluid on scales much less than the Hubble length with negligible 
shear stress, linear perturbations of energy-momentum conservation gives 

1 d f dv\ . 2 
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Here Cs = {dp/ dp)^/"^ is the adiabatic sound speed and the dimensionless entropy 
perturbation is 

cr = — ^ . (3.4) 

p + p 

Equation p.3p says that the density perturbation growth is driven by the Newtonian 
potential $, it is opposed by adiabatic pressure perturbations through the Jeans 
term c^v, and responds to pressure perturbations produced by entropy perturba- 
tions if they are present. It is important to note that equation p.3p is restricted 
to hnear density perturbations. It may be generahzed to incorporate nonhnear per- 
turbations on scales much smaller than the Hubble length, as is commonly done in 
cosmological simulation codes (see Bertschinger 1998 for a review). 

If pressure forces are unimportant on linear scales — as is thought to be valid 
on scales much larger than 10 Mpc, then atomic gas, galaxies, and cold dark mat- 
ter particles all approximately obey equation p.3p with = a = 0. Two key 
assumptions are implicit: 

1. All matter couples universally to gravity; there are no additional long-range 
forces. 

2. Luminous atoms trace dark matter and so do galaxies, or one can correct for 
"biased galaxy formation" (Dekel & Rees 1987; Smith et at. 2007). 

The first assumption involves fundamental physics and is a restatement of the weak 
equivalence principle. It is important to test this assumption on cosmological scales 
but that is beyond the scope of this article (Jain & Khoury 2010). The second 
assumption is astrophysical and represents a serious concern for cosmological tests. 

If these assumptions can be justified, then growth of structure offers an alterna- 
tive means to determining the Newtonian potential A measurement of density 
perturbations ^{a) substituted into equation (13. 3p gives A$. 

(6) Weak gravitational lensing 

Weak gravitational lensing results when the displacements of photon trajectories 
implied by equation (12.61) accumulate over the past light cone so as to deform the 
images of distant sources. The angular defiection arising between and -|- 5zl 
is 

dr 

(5Q; = -Vi($ + *)— (5zL (3.5) 
dzL 

where a is a two-dimensional vector perpendicular to the photon trajectory. The 
gravitational lens equation (Bartelmann & Schneider 2001) 

(3 = 6-— 5a. (3.6) 

relates the observed direction of a ray to the deflection and the direction to 
the source (3 in the absence of deflection. The two distances rs and r^s are the 
comoving angular distances r[x) appearing in the spatial line element "fijdx^dx^ = 

dx^ + r^{x){d9^ + sin^ Od<l)^)- 

The absolute deflection of rays is unobservable since we have no way of knowing 
the true direction of sources in the absence of deflection. However, the gradient 
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of deflection results in image magnification and shear that are observable and are 
described by the inverse magnification matrix 

M- = ^ = 1 + r !^V,V,($ + ^)^dz, . (3.7) 
86 Jo rs dzL 

The two-by-two magnification matrix is parametrized by a convergence k and shear 
parameters 71 and 72: 

M=( ^ + '^ + ^^ 1,^^ ) . (3.8) 

This description is most useful in the weak lensing limit <C 1, 7^ <C 1. Although 
the convergence k cannot be measured directly because the unlcnscd size and flux 
of sources is generally unknown, the shear components 7i can be measured by 
averaging galaxy image distortions over fields of view containing hundreds of distant 
galaxies (Hoekstra & Jain 2008). Thus weak lensing offers a probe of the two- by- two 
traceless Hessian matrix of the sum of the two potentials: 

^ij = (viV,- - i^i, A2) ($ + *) (3.9) 

where A2 = V_l • V_l. Given this matrix as a function of direction, and impos- 
ing appropriate boundary conditions, it is possible to invert the partial derivative 
operator to solve for $ -|- or the convergence k (Kaiser & Squires 1993). 

(c) Power spectra and transfer functions 

An ideal test of gravitation theories would determine ^{x,r) and ^'(03, r) along 
the past light cone giving (with the expansion history) the solution to the grav- 
itational field equations over as much of spacctimc as is accessible to our view. 
Indeed, reconstruction should be the goal. However, many problems impede this ef- 
fort: galaxies (the main cosmological probe) are sparsely distributed through space, 
their properties and small-scale clustering are strongly affected by nonlinear dynam- 
ics and the astrophysics of star formation (e.g., Smith et al. 2009), and building 
up complete samples is an ongoing project for the field of observational cosmol- 
ogy. In the case of CMB anisotropy, accurate measurements arc possible but the 
Sachs- Wolfe and ISW effects cannot be separated from other sources of anisotropy 
generated during recombination and before. Even if they could be separated, the 
ISW effect is most important at low multipolcs where cosmological theories pre- 
dict that the power spectrum has a chi-squared distribution with 21 + 1 degrees of 
freedom (cosmic variance). 

For these reasons, it is useful to measure the statistical properties of the galaxy 
density and peculiar velocity fields, weak lensing shear, and the ISW effect by mea- 
suring correlation functions. Correlation functions can be regarded as a lossy data 
compression technique; they eliminate information about non-Gaussianity but fa- 
cilitate tests of the time- and length-scale dependence of fluctuations and thereby of 
the metric potentials. Each fluctuating field ($, ^, v, u, 7j, etc.) has an autocorre- 
lation with itself and cross-correlations with other fields. The Fourier transform of 
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a spatial correlation function is a power spectrum or cross-power spectrum, which 
itself is related to a correlation in Fourier space. For a Euclidean background, 

^e''= (— ^)p.,(fc) , (3.10) 

where 

{v{k)v{k')) = (2TifP,,{k)6^{k + k') . (3.11) 

These equations assume that the fluctuations are statistically homogeneous and 
isotropic, which holds for quantum fluctuations superposed on a translationally- 
and rotationally-invariant background. 

The time-dependence of the fields and correlations has been suppressed in equa- 
tions p.iop - p.lll) . In linear perturbation theory, all fluctuation fields are given by 
convolutions of a primordial scalar field or possibly a superposition of several such 
convolutions. The simplest category is curvature fluctuations characterized by the 
gauge-invariant spatial curvature perturbation 

n = -^ + Hu, (3.12) 

where u is the total velocity potential (a weighted average over all species) de- 
fined by equation (j3.1bp . The field TZ is constant on scales larger than the Hubble 
and curvature lengths irrespective of the expansion history (Lyth 1985). It is even 
constant on super-horizon scales in modified gravity theories (Bertschinger 2006a) . 

Every other field is given, in Fourier space, by a rotationally invariant transfer 
function multiplying this primary variable: 

v{k,T) = v{k,T)n{k) . (3.13) 

The transfer functions for all matter and metric variables are given by solutions 
of the appropriate field equations with initial conditions TZ{k) — 1. The power 
spectrum of TZ is conventionally made nondimensional by defining 

WMAP 7-year resuhs (Larson et al. 2011) give = (4.93 ± 0.11) x 10"^ for 
k = 0.002 Mpc-i. 

Cross-correlation techniques can be used to extract information about the ISW 
effect that is otherwise hidden by other sources of low-multipole CMB anisotropy. 
Galaxy (mass) density fluctuations at modest redshift correlate with the potential 
fluctuations giving rise to the ISW effect. The cross-correlation of CMB anisotropy 
with galaxy density can therefore extract the ISW signal (Crittenden & Turok 1996; 
Afshordi et al. 2004). 

Cross-correlation techniques can also be used to minimize the systematic effects 
of galaxy bias. Zhang et al. (2007) defined a ratio of weak lensing shear to galaxy 
peculiar velocities based on cross-correlations with the galaxy density. The resulting 
statistic, Eq, can be measured on large scales using weak lensing combined with 
redshift-space distortions in galaxy redshift surveys. 
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4. Potentials in GR and in f{R) theories 

To test a theory of gravity, test the gravitational field equations. As argued in the 
Introduction, the field equations for the unperturbed background are insufficient. 
Therefore we examine the linear perturbations of the field equations. In general 
relativity the field equations are G^^ = SttGT'^j^, with first-order perturbations 
giving 

(A + 3/^)* = aiy , {4.1a) 
i' + n<P = au , (4.1b) 
v]/ _ $ = an , (4-lc) 

-A($ - *) = — 

3 ^ ' p+p 

where 

a = AnGa'^ip + p) . (4.2) 

The right-hand side of equations (j4.ip uses the stress-energy perturbations given in 
equations (13. ip and (|3.2p . An alternative form of equation (j4.1dp is (Bertschinger 
2006a) 



- iCI' ni2^ + ^) + ( 2-H + -H"^ ) $ + i A($ - *) = , (4.1d) 



a d 



cl^^^aS, S = a + ^-^{n\) + ^A7r . (4.3) 



This equation is exact in general relativity with linear perturbation theory for open, 
curved, or flat backgrounds, and for arbitrary matter fields and dynamics. 

If we know nothing about the mass-energy content of the universe, then the 
Einstein equations provide no possible test; we cannot test GR without making 
assumptions about the mass-energy content. We can simply measure the potentials 
and claim there exist matter fields obeying the Einstein field equations. Because 
of the Bianchi identifies, we are even guaranteed that the "stress-energy" tensor 
inferred this way is locally conserved. 

This approach is unreasonable. We aspire to a complete physical characterization 
of mass-energy in the laboratory as well as the cosmos. Although we do not yet know 
the composition of the dark matter and dark energy, the most natural models share 
an important property: the shear stress is small compared with the energy density. 
For example, for linear perturbations of a scalar field or a perfect fluid, tt = 0. 
For free-streaming particles, Att ~ c^z/ where Cs is the sound speed. The largest 
source of shear stress in the standard cosmological model is relativistic neutrinos 
after neutrino decoupling. On scales larger than the Hubble length, and while the 
neutrinos are relativistic (Ma & Bertschinger 1995, with some change of notation), 

2 ip+p)^ , , 

An = —— — — v (4-4) 

5 (p +p)tot 

leading to (Bertschinger & Zukin 2008) 

$ Hp + P)tot ■ ^ ' 
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This ratio is negligible in the matter-dominated era. In order for shear stress to be 
gravitationally important, there must exist a relativistic component which domi- 
nates the mass-energy density of the universe and which has nearly maximal free- 
streaming. Suppose that it is the dark energy. Were such a field present, S in 
equation (I4.3|) would be substantial and would modify the evolution of the grav- 
itational potentials, and thereby would change the evolution of cold dark matter 
perturbations. The observed structure formation could be achieved only if the shear 
stress becomes important after galaxies form (e.g., stressed dark energy, Calabrese 
et al. 2010). Even in this case the ISW effect would be large. 

The combination \[' — $ is called gravitational slip (Daniel et al. 2008; Daniel 
et al. 2009). General relativity with the standard cosmological parameters predicts 
negligible gravitational slip today. 

In alternative theories of gravity the slip is nonzero. For example, in /(i?) theo- 
ries (Sotiriou & Faraoni 2010; De Felice &; Tsujikawa 2010), an additional scalar de- 
gree of freedom is present in the gravitational sector, called the scalaron x (Starobin- 
sky 1980). In equation (|1.3I) . the modification is the term f{R); the scalaron x is 
given by the fluctuations in its derivative. 

The scalaron obeys the wave equation of a field with mass 

The field equations (j4.ip are modified by the addition of terms proportional to x 
and derivatives of fji. For the present purposes the most important modification is 
to the gravitational slip: 

(! + /«)(*-*) -X + a^ ■ (4.8) 

Modified gravity theories, and f{R) in particular, modify the gravitational coupling 
and they introduce a nonzero gravitational slip independently of the matter sector. 

How large is the gravitational slip predicted to be if vr = 0? On scales smaller 
than the Hubble length, a quasi-static solution to the f{R) field equations exists 
with (Hu & Sawicki 2007b; Tsujikawa 2007; Pogosian & Silvestri 2008) 

-I- a^mj 3 

In addition, the Poisson equation (|4.1ap is modified on scales much smaller than 
the Hubble length to become 



$ + ^ 



Combining equations (|4.8p - (|4.9p gives 

$ ~ 4fc2 -I- 3a2m2 



(4.10) 



(4.11) 
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On such scales not only is the gravitational slip nonzero, but the effective gravi- 
tational coupling in the standard Poisson equation A$ = inGeso^Sp is increased 
to 

G f^P + ia^mi 



On scales much larger than the scalaron Compton wavelength m~^, gravity is un- 
modified aside from an overall reduction factor 1 + //j. However, on smaller scales 
the gravitational coupling increases by a factor 4/3 and the gravitational slip be- 
comes ^ - $ = -5$. 

This is a substantial change in the behavior of gravity on small scales, suggesting 
that the f{R) theories might be ruled out based on solar system tests. However, 
the scalaron mass and the coupling factor fn depend on curvature and hence 
on local density. The nonlinearity of the field equations allow for a gravitational 
"chameleon" effect (Khoury & Weltman 2004) that shields gravity in the solar 
system (Hu & Sawicki 2007a) by restricting the rapid change of fa to regions of 
high curvature. It remains to be seen whether star formation and time-dependent 
stellar evolution is compatible with such modifications of gravity in stars. 



5. Current status of tests 

Strong limits on the gravitational slip arc obtained in the solar system using the 
Shapiro time delay (Shapiro 1964), an effect in the propagation of light that is 
sensitive to $ -|- The round-trip time of signals from the Cassini spacecraft near 
Saturn as the signals pass close to the sun is longer than one would expect in flat 
spacetime. The result is an extremely tight limit on the slip (Bertotti et al. 2003): 

- — — — - < 2 X 10 within the solar system . (5.1) 

On galaxy scales, Schwab et al. (2010) have combined stellar dynamics with 
strong gravitational lensing (e.g., the production of Einstein rings) for lenses in the 
Sloan Lens ACS Survey to limit the slip on kiloparsec scales: 

* - $ 

— ^ — — 0.01 ± 0.05 within galaxies on kpc scales . (5-2) 

On scales of a few to a few tens of Mpc, cosmologists have recently been combin- 
ing structure formation measurements with weak lensing (Reyes et al. 2010; Daniel 
et al. 2010; Song et al. 2010) and lensing plus the ISW effect (Bean & Tangmatitham 
2010; Zhao et al. 2010; Daniel & Linder 2010). The Reyes et al. measurement was 
based on the Eq parameter of Zhang et al. (2007), which has the virtue of relative 
insensitivity to galaxy bias and other systematics. Ther results of these tests show 
consistency to date with GR on length scales up to about 50 Mpc. Future surveys 
such as LSST (Abell et al. 2009) and Euclid (Martinelli et al. 2011) will be able to 
improve the limits substantially. 
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